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ABSTRACT 


A  novel  approach  to  battery  lifetime  prediction  has  been 
evaluated  by  application  to  life— cycling  data  collected  for  108 
ESB  EV-106  A-V.  golf  cart  batteries  (tests  conducted  by  TRW  for 
NASA— Lewis).  This  approach  utilized  computerized  pattern 
recognition  methods  to  examine  initial  cycling  measurements  and 
classify  each  battery  into  one  of  two  classes:  "long-lived"  or 
"short-1 i ved".  The  classifier  program  was  based  on  either  a 
linear  discriminant  or  nearest  neighbor  analysis  of  a  training 
set  consisting  of:  each  member  of  the  EV  battery  set  which  had 

failed;  the  relative  lifetime  of  each  member - normalized  with 

respect  to  test  conditions;  and  a  set  of  "features"  based  on 
measurements  of  initial  behavior.  The  raw  data  set  included 
capacity  trends  over  the  first  8  or  9  cycles  and  records  of 
specific  gravity  and  water-added  for  each  cell  after  initial 
cycling.  Features  defined  from  these  raw  data  included  the 
individual  data  items  as  well  as  transf ormat i ons  and  combinations 
of  these  data.  All  features  were  represented  as  standardized 
variables.  It  was  shown  that  lifetime  prediction  of  batteries 
within  the  two  categories  defined  could  be  made  with  about  87% 
accuracy.  It  is  concluded  that  for  a  similarly-manufactured 
battery  set,  relative  lifetime  prediction  could  be  based  on 

initial  measurements  of  the  same  type  examined  here.  _ _ _ _ _ _ — 
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INTRODUCTION 


Traditionally,  battery  lifetime  prediction  has  involved  the 
measurement  of  lifetimes  for  a  sub-set  of  a  "uniform"  population 
of  batteries,  and  then  attributing  the  measured  life 
characteristics  of  the  sub-set  to  the  general  population. 

Under  ideal  circumstances,  one  could  then  predict  the  average 
lifetime,  and  related  variance,  for  a  set  of  batteries,  from  that 
same  population,  operated  under  specified  conditions.  It  has 
not  generally  been  possible,  however,  to  predict  the  lifetime  of 
a  specific  battery,  relative  to  other  batteries  in  the  population. 

For  many  reasons  it  would  be  desirable  to  predict  lifetimes  of 
specific  batteries  in  such  a  way  as  to  discriminate  in  advance 
between  those  which  should  be  "long-lived”  or  "short-lived".  The 
most  obvious  advantage  would  be  to  allow  pre-selection  of  the 
most  reliable  power  sources  for  critical  missions,  such  as  space 
exploration  or  other  remote  operations. 

This  concept  of  specific  lifetime  prediction  was  explored 
previously  by  Byers  and  Perone  (  1  >  for  sealed  Ni /Cd  space  cells 
tested  at  Crane  Naval  Weapons  Support  Center.  The  basic  approach 
involved  the  use  of  pattern  recognition  techniques  to  determine 
if  measurements  of  cell  initial  character i sti cs  could  be  used  to 
predict  the  lifetime  of  specific  cells  relative  to  other  cells 
with  common  origins  operated  under  similar  conditions.  The  basic 
premise  was  that  the  ultimate  fate  of  a  cell  is  reflected  in  a 
multi-variate  examination  of  its  initial  fabrication  and/or 
behavioral  characteristics.  These  measurements  of  initial 
characteristics  become  the  "features"  or  "descriptors”  utilized 
in  pattern  recognition  analysis  to  determine  if  cells  which  are 
destined  to  be  "long-lived"  can  be  discriminated  from 
those  destined  to  be  "short-1 ived". 

The  results  of  this  initial  study  of  pattern  recognition 
lifetime  prediction  (1  )  demonstrated  that  Ni/Cd  cells  from  the 

same  production  lot,  with  similar  fabrication  and  operational 
conditions,  could  be  categorized  from  initial  measurements  with 
virtually  100%  accuracy.  Combinations  of  as  few  as  1  to  3 
features  were  required  to  discriminate  between  predicted 
"short-lived"  and  "long-lived"  cells.  These  features  were  derived 
from  manuf acturer ' s  pre-test  data  documenting  behavior  during 
initial  acceptance  cycle  tests.  The  most  useful  features  involved 
measurements  of  voltage  or  pressure  changes  near  the  end  of  a 
charging  cycle.  Cluster  analysis  of  these  features  suggested  that 
a  quantitative  value  of  relative  lifetime  might  be  assigned  to  a 
specific  cell  based  on  the  average  life  of  its  nearest  neighbors 
in  feature  space.  (This  last  observation  was  very  tentative 
because  of  the  small  size  of  the  individual  clusters  in  the 
limited  data  set.) 
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The  results  o-f  this  initial  study  provide  several  implictions: 

(1)  that  for  a  new  set  of  Ni/Cd  cells  fabricated  identically  to 
the  previous  set  it  would  be  possible  to  predict  lifetimes  of 
specific  cells,  relative  to  all  other  cells  operated  similarly, 
based  on  pattern  recognition  analysis  of  initial  cycle  test  data; 

(2)  that  certain  initial  measurements  may  be  more  sensitive 
lifetime  predictors,  and  these  may  be  useful  in  identifying 
critical  f abrication/operational  factors  dictating  lifetime; 

(3)  that  quantitative  lifetime  prediction  would  be  possible  by 
applying  cluster  analysis  to  a  larger  data  set;  and 

<4)  that  imminent  failure  of  batteries  might  be  predicted  from 
monitoring  of  current  cycle  behavior. 

These  implications  need  to  be  evaluated  in  a  systematic  future 
study.  One  premise  of  the  earlier  work  which  will  be  investigated 
here  is  the  general  applicability  of  the  observations  with  Ni/Cd 
cells  to  other  battery  types.  It  is  the  primary  goal  of  the  study 
reported  here  to  evaluate  this  premise  by  applying  the  same 
lifetime  prediction  techniques  to  a  set  of  Lead/acid  batteries. 

The  rationale  for  examining  the  TRW  data  base  for  life-cycling  of 
108  ESB  EV-106  6-V.  golf  cart  batteries  (2)  is  as  follows: 

the  study  was  wel 1 -designed  and  wel 1 -documented;  the  number 
of  items  with  common  origins  and  test  oonditions  was 
sufficiently  large  for  reliable  pattern  recognition  studies 
(which  require  a  large  ratio  of  patterns  to  features  (  3.  > ) ; 
and  lead/acid  batteries  represent  a  mature  technology  so  that 
positive  results  could  be  directly  useful. 

One  limitation  of  the  use  of  the  TRW  data  base,  however,  is  the 
fact  that  detailed  voltage— time  data  were  not  uniformly 
available.  Thus,  this  study  was  limited  to  an  examination  of  more 
indirect  evidence  of  battery  characteristics.  These  included: 
capacity  trends  over  the  first  8  or  9  cycles;  and  acceptance  test 
measurements  of  specific  gravity  and  water  addition  required  for 
each  cel  1 . 
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Description  o-f  ths  TRW  Data  Base. 

The  test  program  undertaken  by  TRW  (  2  )  Mas  designed  to  apply  a 
daily  charge/di scharge  cycle  program  to  108  lead/acid  6-V. 
batteries  until  failure.  The  conditions  controlled  included 
characteri st i cs  of  a  chopper — controlled  discharge 

(frequency,  duty  cycle);  average/peak  discharge  current;  and  depth 
of  discharge. 

Over  a  2— year  life— cycling  period  69  percent  of  the  batteries 
failed.  Experimental  correlations  showed  that  battery  cycle  life 
was  inversely  proportional  to  depth  of  discharge  and  discharge 
current.  No  significant  effect  on  lifetime  was  detected  for 
different  chopper  discharge  frequencies  and  duty  cycles.  The 
failure  distribution  for  items  with  continuous  (dc)  discharge 
current  was  similar  to  those  items  with  chopper-control  led 
di scharges. 

The  failure  mode  observed  involved  a  gradual  loss  of  capacity  to 
the  half-capacity  failure  point.  23  of  the  failed  batteries  were 
subjected  to  autopsies  which  showed  consistent  evidence  of  cell 
element  aging.  Every  battery  examined  exhibited  short  circuits 
caused  by  metallic  bridging  across  the  plates  at  separator  edges. 
Except  for  two  early  failures,  every  failed  item  examined 
exhibited  buckleo  positive-  pi ates  and  oxidized  positive  grids. 

This  uniformity  of  failure  mode  and  physical  aging  characteristics 
establishes  a  situation  for  pattern  recognition  analysis  which  is 
much  more  nearly  ideal  than  for  the  earlier  study  with  Ni /Cd  cells 
(  1  )  . 


An  Introduction  to  Pattern  Recognition  Methods. 

There  are  many  useful  techniques  for  mathematical  pattern 
recognition.  The  reader  is  referred  to  any  of  several  useful  texts 
on  this  subject  (  4  —  9  )  for  detailed  discussion.  A  brief 
introduction  to  the  concepts  will  be  provided  here. 
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Mathematical  pattern  recognition  methods  take  advantage  o-f  the 
computer’s  ability  to  manage  mul t i -dimensi onal  information  and 
perform  a  series  of  relatively  simple,  but  numerous,  statistical 
and  geometrical  computations.  A  generalized  pattern  recognition 
procedure  involves  several  steps.  The  first  step  involves 
accumulaton  of  observable  data  (d-dimensi onal  pattern  space)  from 
a  physical  system.  Because  the  raw  data  space  may  be  of  large 
dimension,  some  reduction  of  dimensionality  is  desired  to  obtain 
subsequent  reliable  classification.  This  step  involves  the 
definition  of  r — dimensional  feature  space,  where  r  <  d.  The 
reduction  of  dimensionality  should  include  identification  of  those 
features  which  correlate  most  strongly  with  inter-class 
information.  The  next  step  involves  application  of  a  decision 
algorithm  appropriate  for  classifying  the  individual  sources  of 
information  into  any  of  Z  different  classes.  These  decisions  are 
applied  in  r — dimensional  feature  space. 


CLASSIFICATION  METHODS.  Various  general ly-appl icable  mathematical 
procedures  for  pattern  classification  have  been  developed.  Two  of 
these  appear  to  be  particularly  useful  for  the  studies  here.  One 
of  these  involves  Linear  Discriminant  Analysis  (LDA)  <4, 7),  and 
the  other  involves  the  k— Nearest  Neighbor  (kNN)  classification 
criterion  <7, 8). 

Trainable  pattern  classifiers,  a  sub-class  of  learning  machines 
(4),  are  used  in  Linear  Discriminant  Analysis.  The  r  pieces  of 
information  (r  features)  describing  a  pattern  can  be  plotted  as  a 
point  in  f — dimensional  feature  space.  It  is  assumed  that  patterns 
with  similar  properties  will  occupy  the  same  region  of  feature 
space.  LDA  pattern  classification  involves  finding  linear 
boundaries  which  will  discriminate  between  these  spatial  regions. 

A  two-category  pattern  classifier  can  be  defined  by  a  discriminant 
function  which  is  a  scaler,  single-valued  function  of  the  pattern. 
If  the  patterns  to  be  classified  are  linearly  separable,  then  the 
discriminant  function  takes  the  form 

r  +  1 

s  =  .  W|>:j  <1) 

1  =  1 

where  xj  is  the  ith  component  of  a  pattern  having  r  features, 
equals  one,  w*  is  the  weight  corresponding  to  the  ith 
component,  and  s  is  the  scaler  result.  The  category  in  which  a 
given  pattern  is  placed  is  determined  by  the  sign  of  s. 


A  set  of  representative  patterns  of  known  classes,  the  training 
set,  is  presented  sequentially  to  the  classifier.  When  a  pattern 
is  incorrectly  categorized,  the  weights  of  Equation  1  are 
adjusted  in  a  manner  to  correct  the  error.  If  the  training  set  is 
linearly  separable,  this  procedure  will  converge  to  a  single 
weight  vector  which  can  correctly  classify  all  the  patterns. 
Subsequently,  unknown  patterns  of  similar  origins  can  be 
classified  by  the  trained  classifier. 


The  k -Nearest  Neighbor  classification  rule  simply  states  that  an 
unknown  pattern  is  classified  according  to  a  majority  vote  of  its 
k-nearest  neighbors  in  r — dimensional  feature  space. 
Computationally,  the  Euclidian  distances  between  the  single 
r — dimensional  point  representing  the  pattern  in  question  and  all 
other  pattern  points  in  r. -space  must  be  calculated  to  find  the 
nearest  neighbors.  The  distance,  in  r-space,  between  two  points 
i  and  j  is: 


r 


(2) 


Because  the  distances  are  a  non-linear  function  of  the  features, 
the  kNN  method  can  be  applied  to  non-linear  classification 
problems. 


FEATURE  SELECTION  TECHNIQUES.  Feature  selection  involves  reducing 
the  dimensionality  of  a  problem  by  eliminating  pattern  descriptors 
unnecessary  for  classification  and  retaining  a  sub-set  of  pattern 
descriptors  (features)  which  are  required  for  classification. 
Statistical  feature  selection  methods  of  the  principal  components 
type  work  well  for  data  sets  with  well-defined  distribution 
functions. 

For  non-par ametric  classification  problems,  various  transform 
methods  have  been  used  (10  -  13),  as  well  as  a  host  of  other 
basically  empirical  approaches.  However,  most  workers  agree  that 
only  an  empirical  tr i al -and-error ,  al 1 -possible-combi nations 
approach  guarantees  finding  the  optimum  feature  set.  Thus,  in  our 
work  a  systematic  tri al -and-error  feature  elimination  procedure 
was  used  (13),  guided  by  visual  examination  of  feature  plots. 

All -possible-combi nations  of  small  feature  sets  were  also  used. 


RESULTS  AND  DISCUSSION 


Definition  of  Features  and  Sets  for  Lifetime  Prediction. 

LIFETIME  DISTRIBUTIONS.  Figure  1  shows  the  overall  lifetime 
distribution  for  all  10S  batteries  in  the  TRW  study.  (The  large 
block  at  the  upper  end  of  the  distribution  represents  32  batteries 
unfailed  at  the  end  of  the  test  period  of  589  cycles-.)  One  of  the 
most  crucial  steps  in  the  examination  of  the  life-cycling  data  is 
the  assignment  of  battery  lifetime  into  various  categories,  such 
as  "short",  "long",  "average",  or  other.  The  approach  used  in  the 
previous  study  (  1  )  involved  the  use  of  natural ly—occuring 
break(s)  in  the  failure  di stribution (s)  to  define  "long-lived” 
and  "short-lived"  classes.  The  same  approach  was  adopted  here. 

The  validity  of  the  failure  distribution  for  lifetime 
categorization  depends  on  the  uniqueness  of  the  data  sub— set 
selected.  The  set  presented  in  Figure  1  includes  all  items, 
regardless  of  variations  in  test  conditions,  and  as  such  is  not 
very  useful.  In  this  study,  for  example,  it  has  been  observed  that 
lifetime  is  dependent  on  depth  of  discharge  (DOD)  and  on  average 
discharge  current  (IAV).  Figure  2  illustrates  the  differences  in 
lifetime  distributions  for  DOD  —  50  and  DOD  =  75.  (Again,  the 
upper  block  in  each  case  represent  items  unfailed  after  589 
cycles.)  Each  of  these  distributions  could  be  further  sub-divided 
into  sub-sets  with  constant  IAV.  Unfortunately,  these  sub-set 
limits  would  be  too  restrictive  for  analysis  of  the  TRW  data 
base,  because  no  more  than  15  items  could  be  included  in  a 
single  sub-set.  Three  different  nominal  depths  of  discharge  were 

applied - 25,  50,  and  75  percent;  and  7  different  values  of  IAV 

were  employed,  varying  from  20  to  260  A. 

Because  the  reliability  of  pattern  recognition  assignments 
decreases  significantly  when  the  ratio  of  patterns  to  features 
drops  below  about  5  (3),  itis  desirable  to  utilize  a  larger 
pattern  set  to  allow  pattern  recognition  in  a  higher  dimension 
feature  space.  The  basic  approach  taken  here  to  define  larger 
sub-sets  of  the  test  data  for  lifetime  prediction  involved, 
firstly,  normalizing  the  £ycle-life  character i sti c  with  respect  to 
influential  parameters.  Secondly,  the  items  were  grouped 
according  to  DOD  in  examining  the  failure  distributions.  This 
allowed  the  identification  of  3  sub-sets,  correspondi ng  to  the  3 
DOD  values,  25,  50,  and  75.  Because  only  27%  of  those  in  the  first 
sub-set  had  failed  (8  items)  by  the  end  of  the  test  program,  that 
sub-set  was  not  useful  here.  The  other  two  sub-sets  had  81%  and 
97%  failures,  (39  and  29  items  respectively),  and  were  useful 
sub-sets  for  pattern  recognition. 
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NORMALIZED  LIFETIME  DISTRIBUTIONS.  Normalization  of  the 
cycle-life  character i stic  was  accomplished  according  to  one  of 
two  relationships.  The  first  was  referred  to  as  "REGLIFE1*.  It  was 
equated  to  the  deviation  of  observed  lifetime  from  the  regression 
fit  to  the  expression: 

FLCYNO  =  A  +  BilAV  (3) 


where  FLCYNO  is  the  number  of  cycles  to  failure,  and  IAV  is  the 
average  discharge  current.  The  constants,  A  and  B,  are  determined 
from  a  fit  to  all  data  where  DOD  is  constant.  Then,  for  each 
failed  battery  in  the  DOD  sub-set,  the  value  of  REGLIFE  is 
calculated: 


REGLIFE  =  FLCYNO  -  A  -  B*IAV 


<4> 


The  second  normalization  function  was  based  on  the  calculation  of 
a  relative  lifetime,  referred  to  as  "RELIF".  It  was  equated  to 
the  ratio  of  the  observed  lifetime  for  a  specific  battery  to  the 
average  lifetime  for  all  other  batteries  where  DOD  and  IAV  are 
the  same. 

RELIF  =  FLCYNO /AVG (FLCYNO)  (5) 

The  distributions  of  REGLIFE  and  RELIF  were  determined  for  all 
batteries  in  the  two  sub-sets  corresponding  to  DOD  =  50  and  75. 
<For  convenience,  these  will  be  referred  to  as  the  D0D50  and  D0D75 
sub-sets.)  For  the  D0D50  sub-set,  the  maximum  number  of  members 
was  48.  However,  3  batteries  failed  prematurely  at  less  than  173 
cycles  and  were  eliminated  arbitrarily  from  consideration.  This 
provided  a  useful  sub— set  of  45  items.  9  batteries  had  not  failed 
by  the  end  of  the  test  period.  (6  of  these  9  were  tested  under 
the  least  strenuous  conditions  <IAV  =  20  A.)).  Each  of  these 
items  was  assigned  an  arbitrary  lifetime  of  625  cycles,  obtained 
from  an  extrapolation  of  a  plot  of  lifetimes  vs.  IAV.  The  maximum 
measured  lifetime  was  589  cycles,  after  which  the  test  was 
discontinued.  The  median  lifetime  for  the  D0D50  sub-set  of 
45  items  was  495  cycles. 
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For  the  D0D75  sub-set,  the  maximum  number  of  members  was  30. 

There  were  no  premature  failures.  Only  one  battery  was  unfailed 
at  the  end  of  the  test  period.  It  was  assigned  an  extrapolated 
arbitrary  lifetime  of  600  cycles.  The  median  lifetime  for  the 
30- item  D0D75  sub -set  was  403  cycles. 

The  values  of  the  regression  constants  determined  from  the  fits  to 
Equation  3  for  the  D0D50  and  D0D75  sub-sets  are  listed  in  Table  I. 
Figures  3  and  4  present  the  distributions  of  REGLIFE  and  RELIF 
for  the  D0D50  and  D0D75  sub-sets  of  45  and  30  items, 
respectively.  In  contrast  with  Figure  2,  it  appears  in  each  case 
that  the  D0D50  and  D0D75  distributions  might  be  combined  into  a 
single  set.  However,  subtle  differences  in  the  DOD50  and  D0D75 
distributions  are  observed  for  both  REGLIFE  and  RELIF.  One  of 
these  differences  is  obvious  in  comparing  the  ranges  of  the 
distributions  for  the  D0D50  and  D0D75  sub— sets  in  Figures  3  and  4. 

In  both  cases  the  D0D75  range  is  considerably  broader  than  for  the 
D0D50  range.  Another  fundamental  difference  becomes  clear  when  class 
boundaries  are  defined  by  pattern  recognition  as  discussed  below. 


CATEGORIZATION.  Examination  of  Figures  3  and  4  shows  that  some 
natural  breaks  appear  in  the  normalized  failure  distributions. 
These  can  be  used  as  a  first  cut  assignment  of  lifetime  classes 
for  pattern  recognition.  For  example,  in  Figures  3A  and  3B  it 
appears  that  several  passible  "breaks"  might  be  considered  to 
distinguish  between  "short-lived"  and  “long-lived"  classes.  The 
obvious  breaks  in  3A  occur  at  REGLIFE  values  of  ^<-80,  -35,  and 
+70);  those  in  3B  occur  at  REGLIFE  values  of  ''*<-150,  -80,  0,  and 
+100).  However,  several  other  breaks  appear  when  any  regions  of 
the  histograms  are  expanded. 

The  selection  of  any  particular  boundary  for  categorization 
depends  on  two  things.  One  of  these  is  the  purpose  of  the 
categorization.  For  example,  several  different  binary 
categorizations  could  be  considered,  where  the  purposes  might  be 
to  identify  the  very  best  cells;  the  very  worst  cells;  or  to 
simply  distinguish  between  the  two  classes  "better”  and  "worse". 

In  addition,  a  3-class  categorization  could  be  considered  which 
included  the  best,  worst,  and  middle  classes.  Arbitrary 
boundaries  can  be  assigned  based  on  naturally  occuring  breaks 
once  the  purpose  has  been  defined.  For  example,  if  a  simple 
binary  categorization  of  "better /worse"  is  desired,  the  first-cut 
boundaries  for  Figures  3A  and  3B  might  be  -36  and  0,  respectively. 

The  second  criterion  for  selection  of  category  boundaries  is 
based  on  the  observed  performance  of  the  selected  boundary  for 
pattern  recognition  classification  of  cells  with  known 
performance.  The  selected  boundaries  can  be  further  refined  by 
adjusting  them  for  optimum  classification  accuracy  from  pattern 
recognition  examination  of  measurement  features. 
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DEFINITION  OF  FEATURES  FOR  PATTERN  RECOGNITION.  The  -features  used 
•for  pattern  recognition  lifetime  prediction  were  taken  from 
documentation  of  the  preliminary  examination  and  initial 
acceptance  tests  applied  to  all  EV-106  batteries  prior  to  commencing 
life-cycle  testing.  These  included  measurements  of  the  specific 
gravities,  battery  weights,  and  volume  of  water  required  to 
achieve  uniform  levels  for  each  cell /battery  as  received,  as  well 
as  discharge  capacity  values  over  8  or  9  acceptance  cycles. 

These  initial  acceptance  data  were  used  to  generate  pattern 
features  for  each  battery.  The  most  useful  features  fell  into  4 
categories: 

o  Initial  Specific  Gravity 

o  Initial  Water  Volume  Added 

o  Initial  Capacity  Trends 

o  Transformed/Combined  Variables 


A  total  of  lO  features  proved  to  be  useful,  and  these  are 
summarized  in  Table  II. 

A  summary  of  the  variances  of  each  of  the  pattern  features  for 
two  different  sub-sets  is  provided  in  Table  III.  Because  of  the  wide 
disparity  in  values  of  the  features,  all  pattern  recognition  studies 
were  conducted  with  standardized  variables,  where  the  standardized 
•  value,  >:(s),  is  defined: 

x  (s)  =  <x<i)  -  x <ave) ) / (s.d. )  <6> 

Thus,  the  ranges  of  all  standardized  variables  were  ~(+/-3). 

As  expected,  some  of  the  features  were  highly  correlated.  However, 
as  observed  in  previous  studies  (  13  —  15  ),  the  use  of 
statistically  correlated  features  can  be  justified  and. useful  for 
pattern  recognition  where  normal  distributions  are  net  observed. 

This  is  certainly  the  case  for  our  data. 

The  selection  of  an  optimum  feature  set  for  classification 
is  a  crucial  part  of  any  pattern  recognition  study,  and 
the  procedure  used  here  is  described  below. 
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Data  Analysis. 

CLASSIFICATION  PROCEDURES.  Two  different  techniques  were  used  for 
pattern  recognition  lifetime  prediction:  Linear  Discriminant 
Analysis  (LDA) ;  and  k-Nearest  Neighbor  analysis  (kNN).  (These 
were  discussed  in  an  earlier  section.)  The  LDA  method  allows 
accurate  classification  when  classes  can  be  separated  by  a  linear 
boundary  (line,  plane,  hyperplane)  in  feature  space.  Once  a 
discriminant  function  is  found  which  provides  accurate 
classification,  the  application  of  this  function  to  pattern 
recognition  is  computationally  simple.  However,  the  restriction 
to  linearly  separable  classes  precludes  application  in  many  cases. 

The  kNN  method  allows  accurate  classification  as  long  as  items 
of  similar  class  form  clusters  in  feature  space.  These  need  not  be 
linearly  separable  for  classification,  as  long  as  significantly 
different  spatial  distributions  are  obtained.  However,  the  kNN 
algorithm  requires  much  more  extensive  computations  at  the  time 
of  classification. 

For  simplicity,  a  value  of  k  =  1  was  used  for  nearest  neighbor 
calculations.  A  "leave-one-out”  procedure  was  used  to  evaluate 
classification  accuracy.  That  is,  each  item  was  removed  from  the 
training  set,  and  treated  as  an  item  of  unknown  class.  Its  class  is  then 
assigned  to  be  the  same  as  its  nearest  neighbor.  It  is  then  returned  to 
the  training  set,  and  the  next  item  is  removed  for  classification. 

An  iterative  f eature-wei ght i ng  procedure  was  also  used  in  the  kNN 
method  here.  That  is,  feature  weights  were  systematically  varied  by 
multiples  of  2  to  obtain  an  optimum  combination  of  weights  for 
maximum  accuracy. 

FEATURE  PLOTS.  An  examination  of  pattern  distributions  in  feature 
space  provide  useful  insight  to  the  applicability  of  LDA  or  kNN 
classification  techniques.  Figure  5  shows  a  feature  plot  of  INCAP 
and  SGH  for  the  D0D75  sub-set,  where  class  assignments  were  based 
on  the  optimum  boundary  for  the  REGLIFE  distribution  (discussd 
below).  In  this  case  the  two  classes  are  linearly  separated  in 
2-d  feature  space.  Thus,  the  LDA  method  works  very  well  for 
classification.  One  possible  linear  boundary  which  would  provide 
accurate  classification  is  also  illustrated  in  Figure  5.  The  kNN 
method  does  not  work  well  for  this  distribution. 

Figure  6  shows  a  weighted  feature  plot  of  INDL  and  MXH2  for  the 
D0D50  sub-set,  where  class  assignments  were  based  on  the  optimum 
boundary  for  the  REGLIFE  distribution.  Clearly,  the  long-lived 
and  short-lived  classes  exhibit  different  distributions,  but  are 
not  linearly  separable  in  2-space.  When  a  third  feature  (SGH)  is 
added,  the  class  distributions  become  separated  sufficiently  to 
allow  accurate  kNN  classification  (see  results  section). 
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FEATURE  ELIMINATION  PROCEDURF3.  Two  different  methods  were  used  to 
obtain  a  minimum  feature  set  for  optimum  classification.  One  of 
these  involved  using  all  possible  combinations  of  1,  2,  or  3 
features  from  those  defined  in  Table  II.  The  other  involved  a 
sequential  elimination  procedure  (13). 

The  sequential  elimination  procedure  involved  first  conducting 
cl assi f i cati on  with  all  features  used.  Then  one  feature  is  removed 
and  classification  is  again  conducted.  If  classification  accuracy 
is  unchanged  or  improved,  the  feature  is  permanently  eliminated. 

If  not  it  is  returned  to  the  feature  set.  This  process  continues 
until  a  minimal  feature  sub-set  is  obtained  where  no  further 
improvement  in  cl assi f i cati on  accuracy  is  observed  by  elimination. 
This  method  proved  useful  for  minimal  feature  sets  greater  than  3. 
However,  the  elimination  sequence  is  too  arbitrary  to  guarantee 
identification  of  the  optimum  feature  sub-set,  particularly  for 
small  sets.  The  al 1 -possi bl e-combi nat i ons  approach  was  practical 
and  effective  for  up  to  3  features. 


OPTIMIZING  CLASS  BOUNDARIES  FROM  LIFETIME  DISTRIBUTIONS.  Although 
several  different  category  definitions  were  used  in  this  study, 
this  report  focuses  on  the  simple  binary  classification  issue 
where  batteries  were  simply  divided  into  longer-lived  and 
shorter-lived  classes.  It  was  expected  that  this  would  lead  to 
some  overlap  in  the  middle,  and  that  less  than  100'/.  accuracy  would 
be  achieved.  However,  identification  of  these  two  classes  appeared 
to  be  a  realistic  goal  for  real  applications  of  these  methods  to 
battery  lifetime  prediction.  When  only  high  accuracy  was  the  goal 
we  were  able  to  achieve  that  by  discriminating  only  the  very  best 
or  the  very  worst  batteries  from  the  rest. 

The  irethod  used  to  identify  the  optimum  class  boundaries  from 
lifetime  distributions  involved,  first,  selecting 
arbitrarily  a  natural  1 y-occuri ng  break  in  the  distribution. 

Pattern  classification  was  applied  to  the  defined  classes  as  a 
training  set,  and  the  resultant  accuracy  observed.  The  boundary 
was  then  adjusted  in  either  direction  searching  for  maximum 
classification  accuracy. 
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CLASSIFICATION  RESULTS.  Examination  of  Tables  IV  and  V  verifies 
that  pattern  recognition  provides  accurate  prediction  of  lifetime 
class  based  on  battery  acceptance  test  data.  Both  the  LDA  and  kNN 
methods  proved  useful.  The  LDA  method  provided  the  most  accurate 
lifetime  classifications  for  all  sub-sets.  For  each  sub-set, 
overall  classification  accuracy  of  '”85%  was  achieved.  Best 
results  were  obtained  for  the  D0D50  sub-set,  where  classification 
accuracy  for  short-lived  batteries  approached  1007..  High  accuracy 
in  identifying  potential  short-lived  batteries  provides  a 
significant  incentive  for  practical  applications  of  predictive 
lifetime  classification. 

There  did  not  appear  to  be  any  significant  advantage  for  either  of 
the  two  lifetime  normalization  methods,  REGLIFE  or  RELIF.  Both 
worked  well.  However,  the  RELIF  distribution  was  the  only  one  that 
could  be  used  for  the  D0D50/75  sub-set  because  the  REGLIFE 
distributions  for  D0D50  and  D0D75  were  so  different. 

The  size  of  the  D0D50  sub-set  was  varied  in  these  studies  to 
examine  the  effects  of  various  anomolies  in  the  test  items.  These 
are  documented  in  Tables  IV  and  V.  The  primary  concern  was  how  to 
handle  the  batteries  which  had  not  failed  by  the  end  of  the  test. 
In  one  case  (D0D50  sub— set  with  39  items,  (D0D5CX39))  all 
batteries  where  IAV  was  20  A.  were  excluded,  as  most  of  these  had 
not  failed  by  the  end  of  the  test.  The  classification  accuracy 
for  this  sub-set  was  the  highest  of  all,  with  short-lived 
batteries  being  identified  with  1007.  accuracy. 

Another  questionable  test  item  was  a  battery  (s/n  =  16)  in  the 
D0D50  sub-set  which  had  an  exceptionally  low  lifetime  (438 
cycles)  for  an  IAV  of  20  A.  Because  all  other  batteries  tested 
under  these  conditions  were  unfailed  at  the  end  of  the  test  (589 
cycles),  this  battery  might  be  considered  anomolous.  By  way  of 
confirmation,  an  autopsy  (2a)  of  this  battery  revealed  that  the 
negative  plates  were  "hard  and  dry”,  a  condition  not  found  in  any 
of  the  other  autopsies.  Thus,  this  item  was  excluded  from  some  of 
the  sub-sets  examined.  Also,  this  battery  was  excluded  from 
regression  analysis  of  the  D0D50  sub-set.  Thus,  the  REGLIFE 
distribution  is  based  on  the  lifetimes  of  a  44— battery  D0D50 
sub -set . 

The  boundary  values  of  REGLIFE  and  RELIF  required  for  optimum 
classification  accuracy  are  listed  in  Tables  IV  and  V.  It  is 
interesting  to  note  that  for  both  distributions  the  optimum 
boundary  between  short-lived  and  long-lived  cells  is  shifted  to 
larger  values  for  the  D0D75  sub-set  compared  to  the  D0D50 
sub-set.  This  results  in  a  larger  percentage  of  batteries  being 
categorized  as  "short-lived"  in  the  D0D75  sub-set.  This  is  not 
inconsistent  with  the  fact  that  the  actual  depth  of  discharge  was 
'v93%,  as  pointed  out  in  the  TRW/NASA  report  (  2  ).  Moreover,  for 
large  values  of  IAV,  the  effective  depth  of  discharge  approaches 
100%.  Thus,  it  is  not  surprising  that  the  fraction  of  batteries 
which  cluster  together  as  a  short-lived  group  is  larger  for  the 
D0D75  sub-set.  Also,  the  relatively  poor  classification  accuracy 
obtained  when  the  D0D50  and  D0D75  sub-sets  are  combined  is  very 
likely  explained  by  the  significantly  different  distributions  of 
the  two  classes  for  the  two  sub-sets. 
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The  validity  of  the  pattern  recognition  results  is  substantiated 
by  the  low  ratio  of  features  to  patterns  required  for  accurate 
classification  in  each  case.  The  largest  ratio  required  was  for 
the  D0D50<39)  sub-set  with  the  LDA  method  ("“1:8),  (Table  IV). 
Typically,  the  ratio  required  was  'v'l:15.  In  any  case,  the  ratios 
obtained  were  much  lower  than  required  (‘vl:5)  for  credible  pattern 
recognition  classification  (  3  ). 

The  most  useful  features  for  predictive  lifetime  classification 
appeared  to  be  SGH  and  INCAP,  based  on  the  high  frequency  of  their 
appearance  in  the  minimum  feature  sets  for  accurate  classification. 
This  observation  is  certainly  consistent  with  the  intuitive 
perception  that  differences  in  specific  gravity,  water  added,  and 
initial  capacity  trends  should  be  meaningful  predictors  of  battery 
life.  It  is  clear,  however,  that  the  relationships  between 
all  features  studied  and  battery  lifetime  are  non-linear  and 
multivariate. 


CONCLUSIONS 


This  study  has  clearly  demonstrated  the  feasibility  of  predictive 
lifetime  classification  for  uniformly  fabricated  lead-acid 
batteries.  Moreover,  the  utility  of  acceptance  tests  documenting 
trends  in  specific  gravity,  water  added,  and  initial  capacity  has 
been  shown.  The  accuracy  of  predictive  classifications  is 
sufficiently  high,  particularly  for  the  i dent i f i cat i on  of 
short-lived  batteries,  that  the  practical  application  of  this 
method  should  be  explored. 

Perhaps  of  more  importance  is  the  fact  that  this  type  of  study  may 

provide  new  insight  to  factors  which  affect  battery  life - as 

reflected  in  the  useful  features  for  predictive  lifetime 
classification.  For  example,  we  should  like  to  know  why  water -added 
is  a  sensitive  indicator.  To  examine  such  questions,  we  are  currently 
undertaking  a  new  study  documenting  added  water  and  acid  content 
changes  (as  well  as  other  measures)  during  acceptance  cycles  in  a 
life-cycling  experiment  with  lead-acid  batteries. 

The  general  applicability  of  the  predictive  features  and  the 
classification  methods  studied  here  for  lead-acid  batteries  of 
various  origins  remains  to  be  investigated.  In  addition,  it  is 
desirable  to  examine  the  utility  of  more-detailed  charge-discharge 
voltage  trend  data  for  predictive  lifetime  classification. 

Moreover,  it  would  be  desirable  to  examine  a  significantly  larger 
population  of  test  articles  to  evaluate  the  feasibility  of 
quantitative  lifetime  prediction  which  had  been  suggested  in 
the  earlier  study  with  Ni/Cd  cells(  1  ). 
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TABLE  I 

Regression  Constants  -for  Fits  to  Equation  1. 


TEST 

SUB-SET 


REGRESSION  CONSTANTS 
A,  cycles  B,  cycles/A. 


D0D50 


608.9 


-0.7864 


DOD75 


444.7 


-0.2526 


T3.21 
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TABLE  II 

Classification  Features  For  Each  EV-106  Battery 
Based  on  EV-106  Acceptance  Tests 


NAME  TYPE  DESCRIPTION 

SPECIFIC  GRAVITY 

AVSG  AVERAGE  SPECIFIC  GRAVITY  3  CELLS 

WATER  VOLUME  ADDED 

MXH  VOLUME  FOR  CELL  REQUIRING  MOST  WATER 

INITIAL  CAPACITY  TRENDS 

AVERAGE  CAPACITY  OF  ACCEPTANCE  CYCLES 
MAXIMUM  CAPACITY  FROM  ACCEPTANCE  CYCLES 
MINIMUM  CAPACITY  FROM  ACCEPTANCE  CYCLES 

TRANSFORMED  VARIABLES 


AVSG2 

< AVSG) 2 

MXH2 

<MXH) 1 

SGH 

AVSG*MXH 

DLCP 

MXCP  -  MNCP 

INDL 

INCAP/DLCP 

INCAP 

MXCP 

MNCP 
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TABLE  III. 

Feature  Values  for  Two  Sub-sets. 


D0D50 ( 45 )  D0D75  <  30 ) 


FEATURE 

AVG. 

S.O. 

AVG. 

S.D. 

AVSG 

1.278 

0.006 

1.278 

0.004 

AVSG2 

1.633 

0.014 

1 . 634 

0.010 

MXN, ml . 

91.8 

38.7 

91.7 

36 . 4 

MXH2, ml J 

9888 

7728 

9682 

7031 

SGH, ml . 

1 1 . 7x  10* 

4 .  2>:  1 0* 

11.7x10* 

4.6x10 

INCAP, A. h. 

106.9 

1.5 

107.2 

1.  1 

MXCP, A.h. 

114.2 

2.  1 

113.7 

2.2 

MNCP, A.h. 

103.8 

2.  1 

104.3 

1.  1 

DLCP, A.h. 

10.4 

2.3 

9.5 

2.2 

INDL 

10.8 

2.4 

11.9 

2.5 

' 
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TABLE  IV.  SUNNARY  OF  CLASSIFICATION  RESULTS  WITH  LINEAR  DISCRIH INANT  ANALYSIS 
(Long-lived  batteries  ■  dm  1;  short-lived  ■  dm  2) 


DATA  BASE  SUB-SET 


DOO 

NO. 

ITERS 

CLASS 

DISTRIBUTION 

<l)/(2) 

CLASSIFICATION 

CRITERION/BOUNDARY 

FEATURES 

REQUIRED 

CLASSIFICATION 
I  CORRECT  Z(2) 

RESULTS 

1(1) 

SO 

39(a) 

29/10 

RELIF/ (0.909) 

(5)  AVS62,  R3H2,  S6H, 
INDL,  INCAP 

87.2 

100 

82.8 

30 

43 

38/7 

RELIF/ (0.877) 

(4)  HKH2,  S6H,  INDL, 
INCAP 

85.5 

84.2 

85.7 

SO 

44(b) 

32/12 

REBLIFE/ (~34.34) 
A-608.9 

B*-0. 7864 

(3)  S6H,  INDL,  INCAP 

81.8 

92.0 

78.0 

75 

30 

16/14 

RELIF/ (0.977) 

(2)  S6H,  INCAP 

83.3 

85.7 

81.2 

73 

30 

16/14 

RE6LIFE/-15.04) 

(2)  S6H,  INCAP 

83.3 

85.7 

81.2 

50/75 

69(a) 

57/12 

RELIF/ (0.878) 

(4)  N3H2,  INDL, 

DLCP,  INCAP 

81.2 

75.0 

82.5 

(a)  Features  of  batteries  with  IAV  •  20  reaoved  froa  data  bait.  (All  but  one  unfailed  at  end  of  teat.) 

(b)  One  battery  (s/n  «  16)  moved  froa  data  base  because  of  anoaolously  Ion  life  for  IAV  *  20. 


10.13 
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TABLE  V.  SUMARY  OF  CLASSIFICATION  RESULTS  NITH  K -NEAREST  NEIGHBOR  ANALYSIS 
(Long-lived  batteries  *  dau  1;  short-lived  *  dm  2) 


DATA  BASE  SUB-SET 

CLASS 


DOD 

NO. 

ITEMS 

DISTRIBUTION 

(0/(2) 

CLASSIFICATION 

CRITERION/BOUNDARY 

FEATURES 

REQUIRED 

CLASSIFICATION  RESULTS 
l  CORRECT  1(2)  1(1) 

50 

39(a) 

29/10 

REL1F/ (0.909) 

(3)  AVS62,  HXH2,  DLCP 

87.2 

70.0 

93.1 

50 

45 

32/13 

REBLIFE/ (-34.34) 
A»608.9 

B» -0.7864 

(3)  HIK2,  S6H,  INDL 

84.4 

84.6 

84.4 

75 

30 

17/13 

RELIF/ (0.953) 

(2)  SSH,  DLCP  or 

MH2,  DLCP 

66.7 

69.2 

64.7 

75 

30 

18/12 

RE6LIFE/(-22.2> 

A»444.7 

B-0.2526 

(2)  S6H,  INDL  or 

HXH2,  INDL 

76.7 

66.7 

83.3 

50/75 

69(a) 

57/12 

RELIF/ (0.878) 

(4)  AVS62,  HXH2, 

85.5 

50.0 

93.0 

S6H,  INCAP 


(a)  Features  of  batteries  with  IAV  «  20  reeoved  froa  data  base. 
(AU  but  one  (s/n  *  16)  uere  unfailed  at  end  of  test.) 


10.1a 
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Figure  1 


Figure  2 


Figure  3 


Figure  4 


Figure  5 


Figure 


LIST  OF  FIGURES 

.  Lifetime  Distribution  for  All  108  Batteries.  Large  block 
at  upper  end  represents  32  unfailed  batteries  after  589 
cycles. 

.  Lifetime  Distributions  for  the  DOD  =  50  Sub-set  of  48 

items;  and  the  DOD  =  75  Sub-set  of  30  items.  Blocks  at  upper 
ends  represent  unfailed  batteries  after  589  cycles. 

.  Lifetime  Distributions  Based  on  REGLIFE  Normalization 
Function.  A.  For  DOD  =  50  sub-set  of  45  items,  with  3 
premature  failures  excluded.  B.  For  DOD  =  75  sub— set  of 
30  items. 

.  Lifetime  Distributions  Based  on  RELIF  Normalization 

Function.  A.  For  DOD  =  50  sub-set  of  45  items,  with  3 
premature  failures  excluded.  B.  For  DOD  =  75  sub— set  of 
30  items. 

.  Feature  Plot  for  DOD  =  75  Sub— set  of  30  Items.  Lifetime 
class  assignments  based  on  optimum  boundary  from  REGLIFE 
distribution  (see  Table  IV). 

i.  Weighted  Feature  Plot  for  DOD  =  50  Sub-set  of  45  items. 
Lifetime  class  assignments  based  on  optimum  boundary 
from  REGLIFE  distribution  (see  Table  V). 


REGLIFE  DISTRIBUTION  DOD50<45> 


<L  J 

i  i  ■  <■■■  ii  >  JWB 

r^-komTi-roo»j— 


li_ce:luc3zdljjiz:o^— 


i: 


REGLIFE 


REGLIFE  DISTRIBUTION  D0D75<3@) 


llQ£luc3:z>iju:z:o 


REGLIFE 


Figure  4A 


LlC^UJCSIDLlJZEO:*- 


RELIF  DISTRIBUTION  DOD75C30) 


U_Ct:LUC3IDLU2:0^- 


FEATURE  PLOT  DOD75<30) 


CVJ  — «  I  CVI 

I  I 

wZO^ICL 


L D 

o*> 


WEIGHTED  FEATURE  PLOT  00050(45) 


□ 

< BO 


0D 

GO 


LONG  LIFE 
SHORT  LIFE 


472: GAN : 7 1 6 -4 
94 /GEN 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  GEN 


Nd . 

Copies 

Office  of  Naval  Research 

Attn:  Code  413 

800  North  Quincy  Street 

Arlington,  Virginia  22217  2 

ONR  Pasadena  Detachment 
Attn:  Dr.  R.  J.  Marcus 
1030  East  Green  Street 

Pasadena,  California  91106  1 

Commander,  Naval  Air  Systems  Command 
Attn:  Code  310C  (H.  Rosenwasser) 

Department  of  the  Navy 

Washington,  D.C.  20360  1 

Defense  Technical  Information  Center 
Building  5,  Cameron  Station 
Alexandria,  Virginia  22314  12 

r  r.  Fred  Saalfeld 

Chemistry  Division,  Code  6100 

Naval  Research  Laboratory 

Washington,  D.C.  20375  1 

U.S.  Army  Research  Office 
Attn:  CRD-AA-IP 
P.  O.  Box  12211 

Research  Triangle  Park,  N.C.  27709  1 

Mr.  Vincent  Schaper 
DTNSRDC  Code  2803 

Annapolis,  Maryland  21402  1 

Naval  Ocean  Systems  Center 

Attn:  Dr.  S.  Yamamoto 

Marine  Sciences  Division 

San  Diego,  California  91232  1 


No. 

Copies 


Naval  Ocean  Systems  Center 

Attn:  Mr.  Joe  McCartney 

San  Diego,  California  92152  1 

Naval  Weapons  Center 
Attn:  Dr.  A.  B.  Amster, 

Chemistry  Division 

China.  Lake,  California  93555  1 

Naval  Civil  Engineering  Laboratory 

Attn:  Dr.  R.  W.  Drisko 

Port  Hueneme,  California  93401  1 

Dean  William  Tolies 

Naval  Postgraduate  School 

Monterey,  California  93940  1 

Scientific  Advisor 
Commandant  of  the  Marine  Corps 
(Code  RD-1) 

Washington,  D.C.  20380  1 


Naval  Ship  Research  and  Development 
Center 

Attn:  Dr.  G.  Bosmajian,  Applied 
Chemistry  Division 

Annapolis,  Maryland  21401  1 

Mr.  John  Boyle 

Materials  Branch 

Naval  Ship  Engineering  Center 

Philadelphia,  Pennsylvania  19112  1 

Mr.  A.  M.  Anzalone 
Administrative  Librarian 
PLASTEC/ARRADCOM 
Bldg  3401 

Dover,  New  Jersey  07801  1 


472: GAN : 7 1 6 -4 
94/359 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  359 
No. 

Copies 


Dr.  Paul  Delahay 
Department  of  Chemistry 
New  York  University 

New  York,  New  York  10003  1 

Dr.  E.  Yeager 

Department  of  Chemistry 

Case  Western  Reserve  University 

Cleveland,  Ohio  41106  1 

Dr.  D.  N.  Bennion 

Department  of  Chemical  Engineering 

Brigham  Young  University 

Provo,  Utah  84602  1 

Dr.  R.  A.  Marcus 
Department  of  Chemistry 
California  Institute  of  Technology 
Pasadena,  California  91125  1 

Dr.  J.  J.  Auborn 
Bell  Laboratories 

Murray  Hill,  New  Jersey  07974  1 

Dr.  Adam  Heller 
Bell  Laboratories 

Murray  Hill,  New  Jersey  07974  1 

Dr.  T.  Katan 
Lockheed  Missiles  and 
Space  Co. ,  Inc. 

P.  0.  Box  504 

Sunnyvale,  California  94088  1 

Dr.  Joseph  Singer,  Code  302*-! 

NASA-Lewis 

21000  Brookpark  Road 

Cleveland,  Ohio  44135  1 

Dr.  B.  B rummer 
EIC  Incorporated 
55  Chapel  Street 

Newton,  Massachusetts  02158  1 

Library 

P.  R.  Mallory  and  Company,  Inc. 

Northwest  Industrial  Park 

Burlington,  Massachusetts  01803  1 


Dr.  P.  J.  Hendra 
Department  of  Chemistry 
University  of  Southampton 
Southampton  S00  5NH 
United  Kingdom 

Dr.  Sam  Perone 
Chemistry  &  Materials 
Science  Department 
Laurence  Livermore  National  Lab. 
Livermore,  California  94550 

Dr.  Royce  W.  Murray 
Department  of  Chemistry 
University  of  North  Carolina 
Chapel  Hill,  North  Carolina  27514 

Naval  Ocean  Systems  Center 
Attn:  Technical  Library 
San  Diego,  California  92152 

Dr.  C.  E.  Mueller 
The  Electrochemistry  Branch 
Materials  Division,  Research  and 
Technology  Department 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Spring,  Maryland  20910 

Dr.  G.  Goodman 
Johnson  Controls 
5757  North  Green  Bay  Avenue 
Milwaukee,  Wisconsin  53201 

Dr.  J.  Boechler 
Electrochimica  Corporation 
Attn:  Technical  Library 
2485  Charleston  Road 
Mountain  View,  California  94040 

Dr.  P.  P.  Schmidt 
Department  of  Chemistry 
Oakland  University 
Rochester,  Michigan  48063 


No. 

Copies 


1 


1 


1 


1 


1 


l 


1 


l 


472: GAN : 7 1 6 -4 
94/359 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  359 


No . 

Copies 


No. 

Copies 


Dr.  H.  Richtol 
Chemistry  Department 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 

Dr.  A.  B.  Ellis 
Chemistry  Department 
University  of  Wisconsin 
Madison,  Wisconsin  53706 

Dr.  M.  Wrighton 
Chemistry  Department 
Massachusetts  Institute 
of  Technology 

Cambridge,  Massachusetts  02139 
Larry  E.  Plew 

Naval  Weapons  Support  Center 
Code  30736,  Building  2906 
Crane,  Indiana  47522 

S .  Ruby 
DOE  (STOR) 

600  E  Street 

Providence,  Rhode  Island  02192 

Dr.  Aaron  Wold 
Brown  University 
Department  of  Chemistry 
Providence,  Rhode  Island  02192 

Dr.  R.  C.  Chudacek 
McGraw-Edison  Company  . 

Edison  Battery  Division 
Post  Office  Box  28 
Bloomfield,  New  Jersey  07003 

Dr.  A.  J.  Bard 
University  of  Texas 
Department  of  Chemistry 
Austin,  Texas  78712 

Dr.  M.  M.  Nicholson 
Electronics  Research  Center 
Rockwell  International 
3370  Miraloma  Avenue 
Anaheim,  California 


Dr.  R.  P.  Van  Duyne 
Department  of  Chemistry 
Northwestern  University 
Evanston,  Illinois  60201 

Dr.  B.  Stanley  Pons 
Department  of  Chemistry 
University  of  Alberta 
Edmonton,  Alberta 
CANADA  T6G  2G2 

Dr.  Michael  J.  Weaver 
Department  of  Chemistry 
Michigan  State  University 
East  Lansing,  Michigan  48824 

Dr.  R.  David  Rauh 

EIC  Corporation 

55  Chapel  Street 

Newton,  Massachusetts  02158 

Dr.  J.  David  Margerum 
Research  Laboratories  Division 
Hughes  Aircraft  Company 
3011  Malibu  Canyon  Road 
Malibu,  California  90265 

Dr.  Martin  Fleischmann 
Department  of  Chemistry 
University  of  Southampton 
Southampton  509  5NH  England 

Dr.  Janet  Osteryoung 
Department  of  Chemistry 
State  University  of 
New  York  at  Buffalo 
Buffalo,  New  York  14214 

Dr.  R.  A.  Osteryoung 
Department  of  Chemistry 
State  University  of 
New  York  at  Buffalo 
Buffalo,  New  York  14214 


472:GAN:716-4 

94/359 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  359 


Dr.  Donald  W.  Ernst 

No. 

Copies 

Mr.  James  R.  Moden 

Naval  Surface  Weapons  Center 

Naval  Underwater  Systems 

Code  R-33 

Center 

White  Oak  Laboratory 

Code  3632 

Silver  Spring,  Maryland  20910 

1 

Newport,  Rhode  Island  02840 

Dr.  R.  Nowak 

Dr.  Bernard  Spielvogel 

Naval  Research  Laboratory 

U.  S.  Army  Research  Office 

Code  6130 

P.  0.  Box  12211 

Washington,  D.C.  20375 

1 

Research  Triangle  Park,  NC  27709 

Dr.  John  F.  Houlihan 

Dr.  Denton  Elliott 

Shenango  Valley  Campus 

Air  Force  Office  of 

Pennsylvania  State  University 

Scientific  Research 

Sharon,  Pennsylvania  16146 

1 

Bolling  AFB 

Dr.  D.  F.  Shriver 

Department  of  Chemistry 

Washington,  D.C.  203J2 

Dr.  David  Aikens 

Northwestern  University 

Chemistry  Department 

Evanston,  Illinois  60201 

1 

Rensselaer  Polytechnic  Institute 

Dr.  D.  H.  Whitmore 

Troy,  New  York  12181 
< 

Department  of  Materials  Science 

Dr.  A.  P.  B.  Lever 

Northwestern  University 

Chemistry  Department 

Evanston,  Illinois  60201 

1 

York  University 

Dr.  Alan  Bewick 

Downsview,  Ontario  M3J1P3 

Canada 

Department  of  Chemistry 

The  University 

Dr.  Stanislaw  Szpak 

Southampton,  S09  5NH  England 

Naval  Ocean  Systems  Center 

Dr.  A.  Himy 

Code  6343 

San  Diego,  California  95i52 

NAVSEA-5433 

NC  #4 

Dr.  Gregory  Farrington 

2541  Jefferson  Davis  Highway 

Department  of  Materials  Science 

Arlington,  Virginia  20362 

and  Engineering 

Dr.  John  Kincaid 

University  of  Pennsylvania 
Philadelphia,  Pennsylvania  19104 

Department  of  the  Navy 

Strategic  Systems  Project  Office 

Dr.  Bruce  Dunn 

Room  901 

Department  of  Engineering  & 

Washington,  D.C.  20376 

Applied  Science 

University  of  California 

Los  Angeles,  California  90024 

No. 

Copies 


1 


1 


1 


1 


I 


1 


472: GAN: 7 16-4 
94/359 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  359 


No. 

Copies 

M.  L.  Robertson 
Manager,  Electrochemical 
and  Power  Sonices  Division 
Naval  Weapons  Support  Center 
Crane,  Indiana  47522  1 

Dr.  Elton  Cairns 

Energy  &  Environment  Division 

Lawrence  Berkeley  Laboratory 

University  of  California 

Berkeley,  California  94720  1 

Dr.  Micha  Tomkiewicz 
Department  of  Physics 
Brooklyn  College 

Brooklyn,  New  York  11210  1 

Dr.  Lesser  Blum 

Department  of  Physics 

University  of  Puerto  Rico 

Rio  Piedras,  Puerto  Rico  00931  1 

Dr.  Joseph  Gordon,  II 
IBM  Corporation 
K33/281 

5600  Cottle  Road 

San  Jose,  California  95193  1 

Dr.  Robert  Somoano 
Jet  Propulsion  Laboratory 
California  Institute  of  Technology 
Pasadena,  California  91103  1 

Dr.  Johann  A.  Joebstl 

USA  Mobility  Equipment  R&D  Command 

DRDME-EC 

Fort  Belvior,  Virginia  22060  1 

Dr.  Judith  H.  Ambrus 
NASA  Headquarters 
M.S.  RTS-6 

Washington,  D.C.  20546  1 

Dr.  Albert  R.  Landgrebe 
U.S.  Department  of  Energy 
M.S.  6B025  Forrestal  Building 
Washington,  D.C.  20595 


Co 

Dr.  T.  Marks 
Department  of  Chemistry 
Northwestern  University 
Evanston,  Illinois  60201 

Dr.  D.  Cipris 
Allied  Corporation 
P.  0.  Box  3000R 
Morristown,  New  Jersey  07960 

Dr.  M.  Philpot 

IBM  Corporation 

5600  Cottle  Road 

San  Jose,  California  95193 

Dr.  Donald  Sandstrom 
Washington  State  University 
Department  of  Physics 
Pullman,  Washington  99164 

Dr.  Carl  Kannewurf 
Northwestern  University 
Department  of  Electrical  Engineering 
and  Computer  Science 
Evanston,  Illinois  60201 

Dr.  Edward  Fletcher 
University  of  Minnesota 
Department  of  Mechanical  Engineering 
Minneapolis,  Minnesota  55455 

Dr.  John  Fontanella 
U.S.  Naval  Academy 
Department  of  Physics 
Annapolis,  Maryland  21402 

Dr.  Martha  Greenblatt 

Rutgers  University 

Department  of  Chemistry 

New  Brunswick,  New  Jersey  08903 

Dr.  John  Wasslb 

Kings  Mountain  Specialties 

P.  0.  Box  1173 

Kings  Mountain,  North  Carolina  28086 


1 


472: GAN : 7 1 6 -4 
94/359 


TECHNICAL  REPORT  DISTRIBUTION  LIST,  359 

No . 

Copies 


Dr.  Walter  Roth 

Department  of  Physics 

State  University  of  New  York 

Albany,  New  York  12222  1 

Dr.  Thomas  Davis 
National  Bureau  of  Standards 
Polymer  Science  and 
Standards  Division 

Washington,  D.C.  20234  1 

Dr.  Charles  Martin 
Department  of  Chemistry 

Texas  A&M  University  1 

Dr.  Anthony  Sammells 

Institute  of  Gas  Technology 

3424  South  State  Street 

Chicago,  Illinois  60616  1 

Dr.  H.  Tachikawa 

Department  of  Chemistry 

Jackson  State  University 

Jackson,  Mississippi  39217  1 

Dr.  W.  M.  Risen 
Department  of  Chemistry 
Brown  University 
Providence,  Rhode  Island 


No. 

Copies 


Dr.  J.  J.  Brophy 
University  of  Utah 
Department  of  Physics 
Salt  Lake  City,  Utah  84112 


l 


END 

DATE 

filmed 


DTIC 


